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Abstract 

Using CLEO data, we study the production of the anti-deuteron, d, in T(nS) resonance decays 
and the nearby continuum. The branching ratios obtained are ^ dir (T(lS') — > dX) = (3.36 ± 
0.23 ± 0.25) x 10" 5 , B(T(1S) -»• dX) = (2.86 ± 0.19 ± 0.21) x 10~ 5 , and B(T(25) -► dX") = 
(3.37 ± 0.50 ± 0.25) x 10~ 5 , where the "dir" superscript indicates that decays produced via re- 
annihilation of the bb pair to a 7* are removed from both the signal and the normalizing number of 
T(15) decays in order to isolate direct decays of the T(15) to ggg, gg^j. Upper limits at 90% CL 
are given for B(T(4S) -► dX) < 1.3 x 10~ 5 , and continuum production a{e + e~ dX) < 0.031 pb. 
The T(2S) data is also used to extract a limit on Xbj — > dX. The results indicate enhanced deuteron 
production in ggg, gg'y hadronization compared to 7* — ► qq. Baryon number compensation is also 
investigated with the large T(1S) — > dX sample. 
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I. INTRODUCTION 



Anti-deuteron production has been observed in e + e~ collisions at both the T(1S) [l| and 
Z 0] resonances as well as in a variety of other interactions [3[ . The study of anti-deuterons 
rather than deuterons avoids large backgrounds from interactions with beam gas and detector 
material in colliders and nuclear breakup in fixed-target and heavy ion collisions. Since the 
various hadronization processes we wish to explore are expected to be charge-symmetric, 
there is no loss of information incurred by studying only the experimentally cleaner anti- 
deuterons. 

Theoretical descriptions of anti-deuteron formation are generally based on a coalescence 
model, according to which an anti-neutron and anti-proton nearby to each other in phase- 
space bind together [Bl]. Simple calculations may be based on empirical anti-baryon pro- 
duction rates, but subtleties arise. Nearby in phase space largely means nearby in vector 
momentum since the hadronization occurs in a compact region. But the finite size of this 
region and the presence of short-lived intermediate resonances, such as the A quartet, lead 
to questions concerning the necessary degree of coherence, which can only be addressed 
with further assumptions. The combined T(IS'), T(2S) result from ARGUS [l| as well as 
an upper limit from OPAL at the Z ji| were accommodated by Gustafson and Hakkinen 
j^l on the basis of a string model calculation used to supply details of the fragmentation 
process. ALEPH j2| also compared their recent result to this model but limited precision 
and momentum range proscribe any definitive conclusions. A more accurate experimental 
result is desirable to further refine models. 

Practical limitations of particle identification restrict the momentum range over which 
anti-deuterons may be studied. However, the lower mass of the T(IS') means that a larger 
fraction of the momentum spectrum is accessible compared to experiments at the Z pole. 
Also, baryon production in Y(IS') decays is known to be enhanced relative to continuum 
hadronization Q. The Z pole provides a generous rate enhancement but the hadronization 
proceeds via an initial qq pair just as the e + e~ continuum, whereas the T(IS') decays pri- 
marily via three gluons which may be contrasted with nearby continuum qq data. Glue-rich 
T(1S) decays might also produce exotic multi-quark states, beyond qq and qqq [8|. As with 
anti-deuterons, these may form in a similar coalescence process of intermediate hadrons or 
from the primary ggg, gg^ hadronization |9|]. It is also interesting to search for evidence of 
anti-deuteron production inconsistent with coalescence. The frequency with which baryon 
number is compensated via two baryons compared to a deuteron accompanying the anti- 
deuteron may prove useful in this regard. 

Our key result will be a much-improved determination of the rate of anti-deuteron pro- 
duction from Y(IS') — > ggg,ggj hadronization. The momentum dependence of production 
may help discriminate production models and is also used to estimate production outside 
our experimentally accessible momentum range. Given the larger data samples, we do not 
need to combine T(15') and T(2S) production as done |l[ previously, but instead use the 
T(2S) data to limit Xbji^-P) production of anti-deuterons. In addition, we obtain a first 
limit on anti-deuterons from the T(4S*) and an improved limit on continuum production. 

II. DATA AND SELECTION CRITERIA 

We use data collected with the CLEO detector at the Cornell Electron Storage Ring, at 
or near the energies of the T(nS) resonances, where n — 1, 2, 4. The analyzed event samples 
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correspond to integrated luminosities of 1.2 fb 1 on the T(1S), 0.53 fb 1 on the T(2S), 
0.48 fb -1 on the T(4S), and 0.67 fb -1 of continuum data from just below the T(4S). The 
resonance samples contain a total of 21.95 x 10 6 T(1S), 3.66 x 10 6 T(2S), and 0.45 x 10 6 
T(4S f ) decays. 

Smaller effective cross-sections on the T(2S) and Y(3>S) and complications from feed- 
down decrease yields and complicate interpretation of these data. Thus we will emphasize 
the T(IS') sample. The T(2S) sample is used to limit anti-deuteron production from % 6 j(lP) 
decays by assuming that the ggg, ggj production from the T(IS') and T(2S) are identical. 
We choose not to analyze an available T(3S) sample since the statistical error on a branching 
ratio would be quite large. It would also not be possible to separate the contributions from 
Y(IS'), T(25), Xbj(lP) and xw(2P) feed-down from the direct T(35) decays in a meaningful 
way. 

The four innermost portions of the CLEO detector are immersed in a 1.5 T solenoidal field. 
Charged-particle tracking is provided by a four-layer double-sided silicon microstrip detector 
To| and a 47-layer small-cell drift chamber with one outer cathode layer [ll]. The drift 



chamber also provides particle identification via specific ionization (dE/dx) measurements 



Surrounding the drift chamber is a LiF-TEA Ring-Imaging CHerenkov (RICH) detector [12 



followed by a CsI(Tl) calorimeter [13j. Most critical in the current analysis are the drift 
chamber, which covers | cos 6 1 ! < 0.93, and the RICH detector, which covers | cos 6 1 ! < 0.80, 
where 9 is the polar angle with respect to the e + e~ beams. 

Our anti-deuteron track selection proceeds as follows. First, a candidate charged track 
must be consistent with originating from the interaction point. The impact parameter 
with respect to the nominal collision point along the beam direction, Az, must satisfy 
\Az\ < 0.05 m; this distribution is dominated by the physical beam bunch length. The 
impact parameter in the r — cf) plane perpendicular to the beam, Ar, is required to satisfy 
|Ar| < 0.005 m. Since the transverse beam-size is much smaller, the difference here is taken 
with respect to a time-averaged collision point to account for accelerator lattice changes 
and other effects. The collision point can be stable over many days for a fixed lattice. The 
track must be well-measured, based on the reduced \ 2 °f t ne track fit and the fraction 
of traversed drift-chamber layers with good hits. Due to difficulties in reconstructing low- 
momentum tracks (especially considering the large energy loss of the softest anti-deuterons 
in material before the drift chamber) and the shrinking dE/dx separation between anti- 
deuterons and other species at high momentum, we only consider tracks with momenta 
between 0.45 GeV/c and 1.45 GeV/c. We will later estimate the amount of signal outside 
this momentum interval. 

The identification of a quality track as an anti-deuteron relies on the ionization energy loss 
measurement in the drift chamber (dE/dx). To ensure a high-quality dE/dx measurement, 
we only use tracks with at least 10 charge samplings remaining after truncation of the 
highest 20% and lowest 5% of the charge samples for each track. Further, the track angle 
with respect to the beamline, 9, must satisfy | cos^l > 0.2 in order to avoid large gas-gain 
saturation effects present at normal incidence with respect to the chamber wires. This limit 
was chosen by examining the behavior of the large inclusive deuteron sample and observing 
where the success of the corrections applied to compensate for this saturation begin to 
degrade. 

The dE/dx measurement is converted to a normalized deviation 

(dE I ' dx) measure( i (dE I ' dx) expected A /, \ 

Xd = (1) 

° dE/dx 
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with respect to the ionization expected for a real (anti-)deuteron. The dE/dx expected mean 
and resolution (a) include dependencies on velocity (fij = p/m), cos 9, and the number of 
hits used to obtain the measured dE/dx. We accept a track as a deuteron candidate if 
—2 < Xd < +3; the asymmetric cut is chosen to reduce background from the large number 
of protons, which appear at lower values of dE/dx. 

To suppress 7r and p background, we impose requirements on the number of detected 
Cherenkov photons in the RICH detector. Proton suppression is important since they they 
are nearest to deuterons in ionization, while suppression of pions is added since they are so 
numerous; kaon suppression is not employed. For a given particle hypothesis, only photons 
within three standard deviations of the expected ring location are counted; we require fewer 
than five photons for the 7r hypothesis and fewer than three photons for the p hypothesis. 
For the entire momentum range, pions are well above Cherenkov threshold and give more 
than 10 detected photons on average, while protons cross threshold near p = 0.9 GeV/c 
with the mean number of photons increasing with increasing momentum. 

III. YIELD EXTRACTION AND BACKGROUNDS 

Our signal is typified by a well- reconstructed track coming from the interaction point, 
with dE/dx consistent with an anti-deuteron. We choose to use the distribution of the 
normalized deviation, Xd, to determine our raw signal yield. We do this in five 200 MeV/c 
momentum bins spanning 0.45 — 1.45 GeV/c. 

The backgrounds to our anti-deuteron signal are from three main sources. The first is 
particle mis-identification. For most of the momentum range considered, dE/dx separation 
is good; however, since anti-deuterons are very rare compared to the other hadrons, even 
a small resolution or mis-measurement tail may be troublesome. Second, spurious hadrons 
are produced via interactions of beam particles or genuine decay products with residual gas 
in the beampipe or the beampipe and inner detector material. In practice, this is a much 
larger issue for deuterons than anti-deuterons since the gas and material are matter and not 
anti-matter. This is the primary reason we focus on anti-deuterons in this study. Finally, for 
resonance decays, there is a possible non-resonant contribution from the continuum events 
underlying the T resonance peaks. 

Our raw yield and the particle mis-identification background are determined as follows. 
We count the total number of entries between —2 < Xd < +3, denoting this as N. To 
estimate the background from misidentification, we fit the Xd distribution to a Gaussian 
signal peak plus an exponential background shape, as shown in Fig. [TJ The mean and width 
of the Gaussian are fixed from fits to the larger deuteron sample in the data. We then define 
a triangular background shape as shown in Fig. [TJ The apex lies on the background curve 
at the Xd value corresponding to the minimum of the total fit function, between the rapidly 
falling background and the signal peak. The triangle is drawn to decrease to zero height at 
Xd = +3. We denote the area of this triangle between — 2 < Xd < +3 as A. We then take the 
central value of the raw yield as N — A/2 with an error of ±A/4. The rapidly falling fit would 
argue for a lower background, while a small accumulation of events at Xd > +3 balances 
this. In fact, we know little about the background shape other than naively expecting a 
falling shape. Our method spans the range from to A for the background size within ±2(7, 
where a = A/ 4. 

For beam particles or decay products interacting with residual gas or beampipe and 
detector material, the tracks do not peak in both r — and z impact parameters as the 
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FIG. 1: The \d distribution in the momentum range 0.45 — 1.45 GeV/c fit to a Gaussian anti- 
deuteron signal and a falling background function. The background estimation employs the area 
of the triangular region between the dashed cuts, as detailed in the text. 



signal does. We therefore estimate the underlying background from these sources by looking 
at the impact-parameter sidebands. These backgrounds are small and assumed to be flat 
for purposes of the modest extrapolation underneath the peaks. 

The non-resonant contribution is estimated through analysis of our off-resonance data 
sample. Since this contribution is found to be small, we use data taken below the T(4S) 
resonance to subtract it for all T resonances. The off-resonance data is scaled by a factor 
accounting for the ratio of luminosities and the 1/s dependence of the cross section with 
respect to the resonance sample. We will also consider the process T(nS) — > 7* — > qq as 
background, because the physics is the same as the continuum e + e~ — ► qq process, and 
we are interested primarily in the ggg and gg^j decays of the T. Thus, when subtracting 
continuum yields, we adjust the scale factor to account for the additional continuum- like 
events produced via re-annihilation via a 7*. Similarly, the most useful branching ratio will 
be normalized not to the total number of decays, N(T(nS)), but instead to 



N dil = N(T(nS)){l-(3 + R had )B i 



(2) 



counting only the decays proceeding via ggg, ggj hadronization by excluding dilepton de- 
cays, which proceed via a 7* as well. Here, B^ = B(T(nS) — > and Rh a d = 
a(e+e- -> hadrons)/a(e+e-_ i -^ We use R had = 3.56 ± 0.01 ± 0.07 [14] and 
B,.. r 

one. This branching ratio is 



- (2.49 ± 0.02 ± 0.07)% [15J. 
We refer to the resulting branching ratio as the "direct' 
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equivalent to 

B-(T(1S) dX) = S^" 7 ^ - (3) 

£(T(1S) -> ggg, gg-y -> X) 

This will be our central result concerning ggg, gg7 hadronization. For completeness, we also 
present a conventional branching ratio without the modifications described in the preceding 
paragraph. Since we will find that continuum production of anti-deuterons is small, while 
non-ggg, gg7 decays of the T(IS') are significant, the direct branching ratio will be larger 
than the conventional inclusive one. 

Results of the above yield and background determinations are summarized in Table [H 



IV. DETECTION EFFICIENCY 
A. Tracking and RICH Efficiency 

We use Monte-Carlo event samples to study the anti-deuteron efficiency of our tracking 
and RICH criteria. We cannot use anti-deuterons in our simulations since this particle is 
not included in GEANT, which is the basis of CLEO Monte-Carlo software. However, we do 
not expect significant differences between deuteron and anti-deuteron behavior since both 
the RICH detector and the tracking are largely charge- independent, as are our selection 
criteria. Nuclear interactions do distinguish d and d, but given our large final errors, we 
may safely neglect this effect as well. In particular, annihilations in the beampipe or silicon 
vertex detector are estimated to be negligible given our statistics. We also note that our T 
decay hadronization models produce very few deuterons; this leads us to choose the following 
techniques. 

Our first Monte-Carlo sample consists of events with one deuteron and no other detector 
activity. The second consists of overlaying the preceding type of "single-track" events on top 
of a real T(1S) decay from data. The former likely overestimates the efficiency due to the 
quiet detector environment, while the latter likely underestimates it due to excess activity 
(since nothing is removed from the full decay when the signal track is added in). We obtain 
tracking efficiencies of about 70%, with a 10% relative difference between the two methods. 

We average the efficiencies of our two Monte-Carlo samples, taking one-quarter of the 
difference between them as a systematic uncertainty. The resulting efficiency is fairly flat, 
except in the lowest momentum bin of 0.45 — 0.65 GeV/c, where it decreases by about 10% 
of itself. By re-weighting Monte-Carlo events according to the momentum distribution of 
the data across this bin, we find that we are not very sensitive to the detailed spectrum, but 
we do add an additional systematic error for this effect. 

Finally, our signal is consistent with being flat vs. cos8 in the accepted range 0.20 < 
| cos#| < 0.93; we assume it is flat when evaluating the effect of our fiducial cut on the 
track-finding efficiency 



B. dE/dx Efficiency 

The CLEO Monte-Carlo simulation of dE/dx measurements is done at the track-level 
and is based on the calibrated expected means and resolutions. However, (anti-) deuterons 
have not been searched for in any other CLEO analysis to date. Since the calibration is 
quite challenging for the very high ionization of the lower momentum anti-deuterons, the 
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FIG. 2: Deuteron \d distributions (dE/dx normalized deviation) in different deuteron momentum 
bins. 



dE/dx calibration was redone for the data samples used here. These new calibrations offer 
less bias versus parameters such as angle and momentum than the standard versions. But, 
as a result, the CLEO dE/dx simulation designed for the standard calibration is not well- 
suited for our analysis. Instead, we use deuterons from real data, produced by beam-gas 
interactions, to estimate the dE/dx efficiency Our impact parameter cuts ensure that these 
tracks are geometrically similar to signal tracks. 

Fig. [2] shows the deuteron Xd distributions for all five momentum bins; these are mostly 
deuterons from beam-gas interactions or nuclear interactions in the detector. We define 
the dE/dx signal efficiency as tdE/dx = N sig /Nd, where iV S j g is the yield in the interval 
— 2 < Xd < +3 and N d is our estimate of the total number of deuterons for all Xd- We 
estimate N d = (N tot — N tai i/2)±N tai i/4:. Here N tot is the yield in the interval — 5 < Xd < +5 
for the lowest two bins, in — 4 < Xd < +4 for the next two, and in — 3 < Xd < +4 
for the highest momentum bin. N tai i is designed to include a possible tail from the large 
background at low Xd and is taken as the yield in the following momentum-dependent 
intervals: —5 < Xd < — 4 for the first two momentum bins, —4 < Xd < —3 for the next two 
momentum bins, and —3 < Xd < —2 for the last bin. The efficiencies are all about 97%, 
except for the lowest momentum bin, where is it about 88% due to the low-side resolution 
tail. Our systematic error on the dE/dx efficiency comes from propagating the error on N d 
quoted above. 

The width of Xd varies with momentum even after our re-calibration, especially in the 
lowest momentum bin. We determine our sensitivity by re-weighting the momentum distri- 



bution in this bin to better reflect the data, and include an additional systematic error on 
the efficiency. 

C. Systematic Uncertainty Summary 

We now summarize our systematic uncertainties; in each case, we give the range across 
the five momentum bins. The total efficiency uncertainty, including track-finding, selection 
criteria, and yield extraction from the Xd plot, ranges from 4.5% to 16%. In addition to 
systematic issues discussed earlier in this Section, we also considered the agreement between 
data and MC simulations of the number of photons associated with tracks in the RICH 
detector and the stability of results for variations in track-selection criteria. The number 
of T(IS') [T(25)] in our data sample has an uncertainty of 1.4% [1.5%] and the continuum 
luminosity is known to 2%. The resulting total systematic uncertainties range from 6.1% to 
16.0%. These are on average comparable to the statistical uncertainties in the case of the 
T(IS') result, and smaller that statistical errors in all other cases. 

V. RESULTS 

A. Anti-deuteron Production in T(15) 

Table [J shows the observed number of events from T(IS') resonance data, off-resonance 
data, and Ar and Az sidebands in the on-resonance data. After subtracting the latter 
sidebands and properly-scaled continuum contributions, and correcting for efficiency, we 
get the number of d events produced by T(IS') decays shown in Table HT1 The direct yield 
column includes a larger continuum scale factor which accounts for the contribution in which 
bb re-annihilate to a virtual photon and form a qq pair whose fragmentation products contain 
an anti-deuteron. We use this column to get the yield from so-called direct decays mediated 
by ggg and gg'j hadronization. 

TABLE I: Anti-deuteron yields and backgrounds for T(IS') data in momentum bins. 



Momentum (GeV/c) 


On T(15) 


Continuum 


Ar sideband 


Az sideband 


0.45 - 0.65 


60.4 ±7.9 


2.0 ± 1.4 


9.0 ±3.0 


6.0 ± 2.5 


0.65 - 0.85 


77.9 ± 9.2 


1.0 ± 1.0 


2.0 ±1.4 


4.0 ±2.0 


0.85 - 1.05 


71.0 ±8.9 


o o +1 - 2 

u - u -o.o 


1.0 ±1.0 


o o +1 - 2 

u - u -o.o 


1.05 - 1.25 


58.1 ±7.8 


o o +1 - 2 
u - u -o.o 


n n +L2 

u - u -o.o 


o o +1 - 2 
u - u -o.o 


1.25 - 1.45 


46.4 ±7.2 


3.0 ± 1.7 


o o +1 - 2 

u - u -o.o 


2.0 ± 1.4 



To get the anti-deuteron yield in the full momentum range, we fit to the Maxwell distri- 
bution as used in fire-ball models (l6| . 

f(p) = a(3 2 ex V (-E/b), (4) 

where (3 = pc/E, and a and b are free parameters. We include as a systematic uncertainty 
the effect of variations of the shape parameter b within the statistical errors of the fit; we do 
not include any systematic uncertainty for the accuracy of the model itself. The resulting 
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TABLE II: Efficiency-corrected anti-deuteron yields and differential branching rations for T(IS') 
data in momentum bins. 



Momentum (GeV/c) 


Corr'd Yield 


Corr'd Direct Yield 


dB dir /dp (10~ 5 c/GeV) 


0.45 - 0.65 
0.65 - 0.85 
0.85 - 1.05 
1.05 - 1.25 
1.25 - 1.45 


85.6 ± 14.3 
111.3 ±14.2 
106.2 ± 13.8 
92.5 ± 12.6 
60.2 ± 12.7 


82.1 ± 16.4 
109.7 ±15.0 
106.2 ± 14.5 
92.5 ± 13.8 
55.5 ± 14.1 


2.2 ±0.5 ±0.2 
3.0 ±0.4 ±0.2 
2.9 ±0.4 ±0.3 
2.5 ±0.4 ±0.2 
1.5 ±0.4 ±0.2 
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FIG. 3: Momentum dependence of anti-deuteron production in T(IS') decays observed by CLEO 
(filled circles) and ARGUS (open diamonds). The solid line shows the fit described in the text; the 
dashed portion is an extrapolation beyond the momentum range which we observe. 

fit to the CLEO data is shown, along with the earlier ARGUS result, in Fig. [31 Much 
of the CLEO systematic error is correlated point-to-point, but statistics still dominate the 
uncertainty. Note that the ARGUS data extends to higher momentum due to their time- 
of-flight system for particle identification. We also note that ARGUS combined T(1S) and 
T(2S) yields to extract a more precise ggg rate; it is not clear what was assumed about 
possible Xbj production. However, it seems most likely that the difference in our results is 
largely statistical; ARGUS has 19 signal events from both resonances combined. 
The final branching ratio per direct T(IS') — > ggg, gg^y decay is 

B dir (T{lS) -> dX) = (3.36 ± 0.23 ± 0.25) x 10~ 5 . (5) 

For this calculation, we have used only the number of T(1S) which decay via ggg, gg^i as 
our normalization and subtracted a small amount of yield due to bb re-annihilation to 7* 
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based on the observed off-resonance continuum yield. The more inclusive "conventional" 
branching ratio result is 

B(T(15) -> dX) = (2.86 ± 0.19 ± 0.21) x 10~ 5 . (6) 

B. d Production in Y(IS') 

Given that the deuteron signal is expected to be identical to the anti-deuteron signal, 
but with very much larger backgrounds, analyzing for deuterons would not contribute much 
statistically to this analysis. However, we can use deuteron production as a consistency check 
on our anti-deuteron measurement. We make this comparison for the restricted momentum 
range 0.6 — 1.1 GeV/c where the signal-to-noise is best. 

Since none of the backgrounds described above peak in both Ar and Az, we use a sideband 
subtraction to remove them. Empirically, we observe that Ar sidebands are very flat and 
we therefore subtract Ar sidebands from the good Ar sample and fit the resulting Az 
distribution to a Gaussian peak plus a polynomial background. This procedure is displayed 
in Fig. HI Here, we only use deuterons which satisfy —2 < Xd < +3 and ignore the small 
backgrounds from other particle types. The resulting deuteron yield is 352.8 ± 88.6, about 
1.7 standard deviations from the anti-deuteron yield of 201.0 ± 14.2 in the same 0.6 - 1.1 
GeV/c momentum range. The Az width is somewhat narrower but similar to that for 
anti-deuterons. 

C. Discussion of d Baryon Number Compensation 

Another way to explore consistency of the d and d yields is to employ baryon number 
conservation. Assuming many of the d (d) events are compensated by pp or pn (pp or pn), 
requiring at least one p (p) in the event may decrease background appreciably. As shown 
in Fig. [5] and summarized in Table IHH after the standard selection criteria, we begin with 
13140 deuterons and 338 anti-deuterons candidates (signal plus background). Our proton 
identification requirements are: \x P \ < 4 for 0.30 — 0.85 GeV/c, and \x P \ < 3 for 0.85 — 1.15 
GeV/c. Here, Xp is a normalized dE/dx deviation with respect to the proton hypothesis 
and we do not accept candidates outside the two contiguous momentum windows indicated 
above. With this definition of protons, we can study the effect of cuts on the number, n p , of 
protons (anti-protons) in anti-deuteron (deuteron) events. If we require n p > 0, we are left 
with 149 d and 898 d; while the non-decay deuterons are very much reduced, the asymmetry 
indicates residual background from random anti-protons not associated with any baryon 
number compensation. This is further verified by examining the Ar, Az distributions in the 
fourth row of Fig. [51 The excess is consistent with a spurious deuteron in coincidence with a 
real physics event containing the anti-proton. A fit to the Ar distribution yields 122.8 ±16.9 
events for the sharp peak, which is now consistent with the anti-deuteron yield. Adding a 
requirement of n p > 2, we obtain 31 d and 35 d, which are quite consistent with equality, 
implying that most spurious d have been removed. (Note that we never observed n p > 2). 
The remaining d events peak well both in Ar and Az, as expected. 

Returning to just the clean anti-deuteron sample, we can take a more quantitative look 
at baryon number conservation. By considering the effect of proton-finding efficiency we can 
determine approximately how baryon compensation is distributed among pp, pn, np, nn; we 
will consider compensation by a d later. 
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FIG. 4: Anti-deuteron and deuteron sample distributions for the momentum range 0.6 — 1.1 GeV/c 
in the T(IS') data. Top row: Ar and Az for anti-deuterons. Middle row: Ar and Az for deuterons. 
Bottom left: Az for deuterons in Ar sidebands. Bottom right: Fit to Az for deuterons after 
subtraction of the Ar sidebands in the previous panel. The signal Gaussian width is fixed to the 
anti-deuteron signal width from the fit to data in the upper right panel. 



TABLE III: Anti-deuteron and deuteron yields from T(15) data with requirements on accompa- 
nying protons and anti-protons. 



Standard Cuts, plus: 


# anti-deuteron candidates 


# deuteron candidates 




338 


13140 


> p 1 > Op 


149 


898 


> 2p 1 > 2p 


31 


35 



In the —2 < Xd < +3 anti-deuteron signal region, we observe 338 candidate events from 
the T(1S) data sample. Each of these events contains only 1 anti-deuteron. Among these 
338 events: 189 events contain no protons, 118 events contain 1 proton, and 31 events 
contain 2 protons. Assuming d is compensated by pp, pn, np, nn (neglecting d for now) 
with an equal probability of 25%, we may estimate what is expected, given a proton-finding 
efficiency. We cannot distinguish pn and np, but it makes the assumed equality clearer to 
list them separately. For the proton identification cuts given above, the efficiency is about 
60%, where we assume the spectrum of protons accompanying deuterons is similar to the 
inclusive proton spectrum. Part of the loss of efficiency is due to protons with momenta 
outside our accepted range of 0.30 — 1.15 GeV/c. 
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FIG. 5: d and d yields in the momentum range 0.45-1.45 GeV/c with additional criteria as de- 
scribed. Top row: \d for deuterons (left) and anti-deuterons (right) with all standard cuts, includ- 
ing —2 < Xd < +3. Second row: As above, but requiring at least one anti-proton (left) or proton 
(right) candidate in addition. Third row: As above, but requiring two anti-proton (left) or two 
proton (right) candidates instead. Fourth row: Ar and Az for deuterons with n p > 0. The curve 
on the left shows a fit to a signal Gaussian plus a polynomial background. Fifth row: Ar and Az 
for deuterons with n v = 2. 



Folding in this approximate efficiency, we predict 30 events containing 2 protons, with 31 
observed, and 142 events containing 1 proton, with 118 observed, as summarized in Table HVl 
Within the limits of our uncertainties and assumptions, our data is consistent with baryon 
number conservation occurring with roughly equal probabilities for accompanying pp, pn, 
np, or nn. 



TABLE IV: Anti-deuteron yields from T(15) decay for various numbers of accompanying protons. 





All 


#p = 


#p = l 


#p = 2 


Observed 


338 


189 


118 


31 


Predicted 




166 


142 


30 



It is also interesting to look for compensation of an anti-deuteron by a deuteron; this is 
found to occur at the 1% level. Fig. [H] shows one of our four possible ddX events, which 
is nearly fully reconstructed. Inspection of these four dd candidate events reveals that one 
of them is consistent with a through-going deuteron track (presumably from a cosmic-ray 
interaction) faking a dd pair. The remaining three are consistent with true dd. Through- 
going deuterons might constitute a non-negligible background to our anti-deuteron yield, if 
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FIG. 6: An T(15) — > 3tt + 3ir~ ddX candidate event in the CLEO detector, viewed along the beam 
axis. The d track (p = 0.55 GeV/c) is highlighted. The d track (p = 0.84 GeV/c) is the one 
with most energetic calorimeter shower; the circle size is proportional to the shower energy. The 
difference between the center-of-mass energy and the total energy of observed particles, based on 
tracking and particle identification, is about 120 MeV. 

the inward deuteron track passed our anti-deuteron cuts, but the outgoing deuteron fails. 
Therefore, we have searched for such events with relaxed cuts. We find none in the Ar 
and Az sidebands, nor do we find any anti-deuteron candidate events where there is a 
lower-quality track candidate failing our cuts back-to-back with our candidate track. We 
conclude that this faking mechanism is rare and the one event seen was a somewhat unlikely 
occurrence for our data sample size. 

D. d Production in 2S, 4S and Continuum 

We now summarize results from other T resonances and the continuum. 

In T(2S) data, 69 anti-deuteron events are observed. This sample has the same back- 
ground sources as the T(IS') data, but contains several possible sources of anti-deuteron 
signal. These include: (i) T{2S) -> T(IS) X, followed by T(lS) -> dX', (ii) T(2S) -> 
999-1 991-1 an d (hi) T(2S f ) — > jXbjX, followed by XbJ dX'. We may subtract process (i) 
based on known branching ratios. Separating (ii) and (iii), for example by looking for the 
transition 7 in (iii), is not feasible with our limited statistics. However, we can assume that 
the rate for the direct decay (ii) is equal to the analogous T(IS') process, and look for any 
excess from Xcj decays. This is interesting since the Xcj decay via gg for J = 0, 2 and via 
gqq for J = 1 and thus access distinct hadronization processes. 
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After background subtractions analogous to the T(1S) case, we find 58.3 ± 8.6 signal 
events, which translates to 

B(T(2S) ^d + X) = (3.37 ± 0.50 ± 0.25) x 10~ 5 . (7) 

To isolate this rate, we subtract contributions from the processes e + e~ — > T(2S) — > 7tttT(IS) 
and e + e _ — > T(25) — > jjT(1S) (two-photon transitions via the Xbj states) assuming that 
these processes dominate inclusive T(1S) production. We must further assume that direct 
ggg, gg^ decays of the Y(2S) produce antideuterons at the same rate as the T(IS'). We are 
left with an insignificant excess, and extract a 90% CL upper limit for a weighted average 
of the Xbj states of 

£(£(T(2S) - 7Xw(LP)) x BilxUlP) -> dX))/J2B(T(2S) - 7 *6 j (LP)) < LI x 10" 4 . 
j j 

(8) 

This limit is not stringent enough to draw firm conclusions on anti-deuteron production in 
these distinct gg and gqq hadronization processes in contrast to ggg, ggj. 

In T(4S) data, 3 d candidates are observed. Based on Ar and Az sidebands and the 
continuum data, we expect 5.2 background events. For both this limit and the following 
continuum production limit, we ignore any possible backgrounds in the Xd distribution. We 



obtain a 90% CL upper limit, using the Feldman- Cousins method [17J, of 



B{T(AS) -> dX) < 1.3 x 10" 5 . (9) 

This limit is not very stringent in view of the dominance of BB decays of the T(45). 

A 90% CL upper limit result for continuum production is also obtained, based on 6 events 
with 1.5 expected background: 

a{e + e~ -> dX) < 0.031 pb, at ^ = 10.5 GeV. (10) 

Given that the continuum hadronic cross-section at \fs = 10.5 GeV exceeds 3000 pb, we see 
that fewer than 1 in 10 5 qq hadronizations results in anti-deuteron production, noticeably 
less than for ggg, gg'j hadronization. 



VI. CONCLUSIONS 

Using CLEO data, we have studied anti-deuteron production from T(nS) resonance de- 
cays and the nearby continuum. T(IS') and T(2S) production rates are presented separately 
for the first time and combined to limit anti-deuteron production from ^j(IP) states. First 
limits on production from the T(45*) and an improved continuum production limit are given. 

The results confirm a small but significant rate from hadronization of T(nS) — > ggg, gg^ 
decays, for n = 1,2. However, no significant production from qq hadronization is observed; 
our qq limit is more than three times smaller than the observed rate from ggg hadronization. 
Thus, the results indicate that anti-deuteron production is enhanced in ggg, gg'j hadroniza- 
tion relative to qq. 

We observe that baryon number conservation is accomplished with approximately equal 
amounts of accompanying pp, pn, np, nn. We also found three dd events; it is not imme- 
diately clear if double coalescence from initial baryons and anti-baryons can accommodate 
this rate, or if this is evidence for a more primary sort of (anti-)deuteron production in the 
hadronization process. 
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